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ABSTRACT 

Making robust predictions for the phase space distribution of dark matter at the solar 
neighbourhood is vital for dark matter direct detection experiments. To date, almost 
all such predictions have been based on simulations that model the dark matter alone. 
Here, we use three cosmological hydrodynamic simulations of bright, disc dominated 
galaxies to include the effects of baryonic matter self-consistently for the first time. 
We find that the addition of baryonic physics drastically alters the dark matter profile 
in the vicinity of the Solar neighbourhood. A stellar/gas disc, already in place at high 
redshift, causes merging satellites to be dragged preferentially towards the disc plane 
where they are torn apart by tides. This results in an accreted dark matter disc that 
contributes ^ 0.25 — 1.5 times the non-rotating halo density at the solar position. The 
dark disc, unlike dark matter streams, is an equilibrium structure that must exist in 
disc galaxies that form in a hierarchical cosmology. Its low rotation lag with respect to 
the Earth significantly boosts WIMP capture in the Earth and Sun, boosts the annual 
modulation signal, and leads to distinct variations in the flux as a function of recoil 
energy that allow the WIMP mass to be determined. 
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1 INTRODUCTION 

The case for dark matter in the Universe is based on a wide 
range of observational data, from galaxy rotation curves 
and gravitational lens ing, to the Cosmic Microwave Back- 
ground Radiat ion (e.g. dc Blok ct al. 2001': (Read et al.ll2007l : 
iDunklev et "all [2008; Clowe ct al 2006). Of the many plau- 
sible dark matter candidates in extensions to the Stan- 
dard Model, Weakly Interacting Massive Pa rticles (WIMPs) 
stand out as well-motivated and detectable l|jungman et al.l 
Il996t ). giving rise t o many experiments desig ned to detect 
WIMPs in the lab llBernabei et al.,.200a: Aprile et al.ll2005l : 



IVogelsberger et al] |2008| ). These exquisitely describe the 
amount and properties of dark matter substructure. But 
they do not model the baryonic component of the Milky 
Way that pr esently dominates the mass interior to the 
solar radius (iDehnen fc BinnevI Il998l : iKlvpin et all |2002| : 
IWidrow fc Dubinskill2005i r and likely did so since redshift 
z = 1, when the mean merger rate in a A CDM cosmology 
peaked ([Da lcanton fc Bernstein 2002; Bc nsbv et al.l l2007l : 
iGottlober et al.ll200ll : Ikampczvk et al.ll2007f r 



ICDMS Collaborationll2008l '). Predicting the flux of dark mat- 
ter particles through the Earth is key to the success of such 
experiments, both to motivate det ector design, and for the 
interpretation of any future signal l|jungman et al.l|l996l ). 

Most previous work has modelled the distribution of 
dark matter at the solar neighbourhood using, or extrap- 
olating from, cosmological simulations that model only 
the dark matter (Stadcl ct al. 2008 
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In recent work, we demonstrated that the Milky Way 
stellar/gas disc is important because it biases the accretion 
of sate llites, causing th em to be dragged towards the disc 
plane (|Read et al.l 120081 '). As these satellites are torn apart 
by tidal forces, their accrete d material settles into a thic k 
disc of stars and dark matter (|Read et al. l l2008l : lLake|[l989l ). 
The dark disc, unlike dark matter streams, is an equilibrium 
structure that must exist in disc galaxies that form in a 
hierarchical cosmology. Its low rotation lag with respect to 
the Eart h significantly boosts WIMP capture i n the Earth 
and Sun (|Bruch et al.i2008d : lB"ruch et aLlbOOSal ). boosts the 
annual modulation signal, and leads to distinct variations in 
the flux as a function of recoil energy that a llow the WIMP 
mass to be determined (|Bruch et al.ll2008bl ). 
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In [Read et al.ll|2008l ). we used dark matter only simula- 
tions to quantify the expected merger history for Milky Way 
mass galaxies in ACDM, and a series of isolated coUisionless 
satellite-Milky Way merger simulations to estimate the ex- 
pected properties of the dark disc. We found that, for plau- 
sible merger histories, the dark disc contributes ~ 0.25 — 1 
times the halo density at the solar neighbourhood. In this 
paper, we make the first attempt to include the baryonic 
matter self-consistently, using three ACDM cosmological hy- 
drodynamic simulations of Milky Way mass galaxies. Both 
approaches are complementary in quantifying the expected 
properties of the dark disc. Our previous approach allowed 
us to specify precisely the properties of the Milky Way disc 
and its merging satellites at high redshift; our current ap- 
proach is fully self-consistent. We model the radiative cool- 
ing, star formation and feedback physics that lead to the 
formation of realistic disc galaxies in their cosmological con- 
text. Recent improvements in such simulations have made 
it possible to form galaxies with a significant disc com- 
ponent, and reprodu ce some basic galaxy properties over 
a range of rcdshifts (" Zavala et al.' 2008; 'Scannapioco et al.l 
[2OO8 ; Govcrnato ct al. 2007; Pontzcn et al. 2008). However, 
due to the stochastic nature of galaxy assembly in ACDM, 
and our limited sample of simulated galaxies, we do not at- 
tempt to form genuine Milky Way analogues. Instead, we 
choose three galaxies with mass comparable to that of the 
Milky Way that span a range of interesting assembly histo- 
ries. We infer from these plausible dark disc properties for 
the Milky Way. 

This paper is organised as follows. In 321 we describe the 
state-of-the-art cosmological hydrodynamic simulations. In 
^ we present our results. In ^ we discuss our numerical 
limitations and the implications of our results for the Milky 
Way. Finally, in ^ we present our conclusions. 



2 THE COSMOLOGICAL HYDRODYNAMIC 
SIMULATIONS 

We use three cosmological hydrodynamic simulations of 
Milky Way mass galaxie s, two of which (MW1,H258) have 
already been presen ted in'B rooks et al.l ll2007l).lBrooks et al.l 



20081), iMaver et al. (2008). iGovernato et alj (120081 ), aiid 



iGovernato et al.' (2008). A ll thre e were run with the GASO- 
LINE codeHWadsley et al.ll2004l ) using "blastwave feedback" 
(|Stinson et al. 20061 ): the simulation labels, parameters and 
choice of cosmology are given in Table[T] We also ran a fourth 
simulation, H258dark. This had the same initial conditions 
as H258, but was run with only dark matter particles, and at 
slightly lower mass resolution. (Note that since the matter 
density of the Universe Q.m is fixed, the dark matter density 
is higher in H258dark than in II258.) The final outputs were 
mass and momentu m centred us i ng the 'shrinking sphere' 
method described in iRead et al.l (|2006al ) , and rotated into 
their moment of inertia eigenframe with the 2 axis perpen- 
dicular to the disc. For MWl, H204 and II258, the eigen- 
frame was found using the stars, for H258dark it was found 
using the dark matter. 

The three galaxies MWl, H204 and H258 were chosen 
to span a range of interesting merger histories, rather than 
as perfect Milky Way analogues. MWl had a very quiescent 
merger history, with no major mergers after redshift z = 2; 



H204 had several massive mergers after redshift z = 1; and 
H258 was an extreme case with a very massive ^1:1 merger 
at 2 = 1. In 33 we discuss which of these is most like our 

own Galaxy. 

The analysis was performed as in iRead et all (|2008h . 
The subhalos inside each 'Milky Way' and at each redshift 
output wer e identified usin g the Amiga Halo Finder (AHfQ) 
algorithm IG iW et aklliooj ). The subhalos were then traced 
backwards in time, the progenitor of each being the subhalo 
at the previous redshift that contains the majority of its par- 
ticles. To avoid ambiguities in this halo tracking, we ordered 
subhalos by mass so that subhalos were linked to their most 
massive progenitor not already claimed by a larger subhalo. 
A final complication can occur if two subhalos are about to 
merge. In this situation AHF sometimes over-estimates the 
mass of the smaller of the two. We dealt with this problem 
by searching for sudden spikes in mass at pericentre and 
removing these by assigning instead the mass found at the 
previous output time. We define a subhalo as merged (dis- 
rupted) if it has less than a tenth of its peak circular speed 
considered over all times; our results are not sensitive to this 
parameter. 



3 RESULTS 

3.1 Forming a disc of dark matter 

Figure [lla-c) shows the distribution of rotational velocities 
at the solar neighbourhood {7 < R < 8kpc; \z\ < 2.1 kpc) 
for the stars (red) and dark matter (black) in MWl, H204 
and H258. Overplotted in panel (c) is the distribution for 
H258dark, modelled with dark matter alone (black dotted). 
Our results are not sensitive to the position of the 'solar 
neighbourhood' slice (a choice of 5.5 < R < 6.5 kpc gave 
near-identical results). With the exception of H258dark, all 
of the dark matter distributions are skewed towards the 
stars. Figure [2ja-d) shows a decomposition of these dis- 
tributions. All three baryonic simulations require a double 
Gaussian fit to their dark matter distribution (smooth 
black and black dotted lines). The second Gaussian rotates 
rapidly and matches well the material accreted from the sum 
of the four most massive satellites (black dotted histogram), 
though our fit did not require this. (Nor did our fit require 
either Gaussian to have mean rotational velocity = 0.) It 
is this accreted, rapidly rotating material that we call the 
'dark disc'. 

Figure [2ja-d) shows how the material accreted from the 
four most massive merging satellite^ (green, blue, magenta 
and cyan) contributes to the dark disc for MWl, H204 and 
H258. Figure [21 e-h) shows the decay in radius as a function 
of redshift z of these satelUtes; and Figure Eli-l) shows the 
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Simulation 




{Na^,N„Ngas)/lO'^ 


min(Afd„, A^., Mgas)/105 Mq 


^ dm ,*,gas / C 


^,<300kpc ,-^Ql2 

dm 1 ^ 


MWl 


(0.3,0.7,0.9,0.7) 


(2.8,3.1,1.5) 


(7.6,0.2,0.3) 


0.3 


1.1 


H204 


(0.24,0.76,0.77,0.73) 


(4,3.3, 1.7) 


(10.1,0.41,0.58) 


0.35 


0.8 


H258 


(0.24,0.76,0.77,0.73) 


(3.5,2.2,1.4) 


(10.1,0.35,0.58) 


0.35 


0.75 


H258dark 


(0.24,0.76,0.77,0.73) 


(3.5,-,-) 


(12.25,-,-) 


0.35 


0.9 



Table 1. Simulation labels and parameters. From left to right the columns show the simulation label, the cosmological parameters used, 
the number of dark, star and gas particles at redshift 2 = 0, the minimum dark matter, star and gas particle masses at 2 = 0, the dark 
matter, star and gas force softenings (these are always equal), and the dark matter mass within 300 kpc at 2 = 0. H258dark was set up 
with the same initial conditions as H258, but run with only dark matter particles, and at slightly lower mass resolution. 
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Figure 1. (a-c) The distribution of rotational velocities at the solar neighbourhood (7 < R < 8kpc; \z\ < 2.1 kpc) for three simulated 
Milky Way mass galaxies MWl, H204 and H258. The lines show the dark matter (black) and stars (red). The dark matter distribution 
for H258dark, simulated with dark matter alone, is overplotted on (c) (black dotted). 



decay in angle 9 to the host galaxies' disc as a function of 
redshift 2. 

The mass and rotation speed of the dark disc increase in 
the simulations with more late mergers. MWl has no signif- 
icant mergers after redshift z — 2 and has a less significant 
dark disc, with rotation lag with respect to the stars (red 
lines in Figure[T|) of ~ 150km/s, and dark disc to halo den- 
sity ratio of Pddisc/phalo = 0.23 (obtained from the dou- 
ble Gaussian fit). H204 and H258 both have extreme dark 
discs with Pddisc/phalo > 1 and rotation lag with respect 
to the stars of ^ 60km/s; they both have massive mergers 
at redshift z < 1. 

Figure [2{i-j) demonstrates that disc plane dragging is 
responsible for the formation of the dark disc. In MWl, the 
green satellite is dragged towards the disc plane, the ma- 
genta and cyan satellites start out close to the disc plane, 
and the blue satellite merges at high inclination angle. Fig- 
ure [2Ja) shows the contribution to the dark disc owing to 
each of these satellites. The magenta satellite contributes 
the most, being both low-inclination and massive, then the 
green. The cyan satellite is of too low mass to contribute sig- 
nificantly, while the blue satellite contributes little rotating 
material because of its high inclination. These results con- 
firm our expectat ions from isolated disc-satellite merger sim- 
ulations (jRead et al.. 2008). Similar results can be seen for 
the four most massive mergers in H204. Although initially 
on high inclination orbits, the magenta and cyan satellites 
complete enough peri-centre passages to be dragged down 
into the disc plane and contribute significantly to the dark 



disc; the blue and green satellites also contribute in equal 
measure, though somewhat less than the magenta and cyan 
satellites owing to their higher final inclinations. H258 ap- 
pears to present a similar picture. The green satellite is a 
near merger that starts out near the disc plane and 

contributes nearly all of the dark disc. However, mergers 
of this mass ratio define the post-merger plane of the disc. 
They will lead to highly rotating, albeit hot, dark matter 
discs, whatever their initial inclination. Yet there is no dark 
disc in H258dark that has the same merger, but no 

baryonic material. This suggests that a second mechanism, 
extra to disc-plane dragging, is important for the formation 
of dark discs. We discuss this below in H3.2I 

It is interesting that in all three galaxies, none of 
the four most massive satellites contribute significant ret- 
rograde material. We will investigate this further in fu- 
ture work, but note here that retrograde mergers are 
suppresse d both because of reduced dyna mical friction 
( Quinn fc GoodnianI 19861 : Read et al] |2008| ). and reduced 
tidal forces (|Read et al.ll2006bl ). 



3.2 



Maintaining the dark disc: the importance of 
halo shape 



Figure [3je-h) shows projected density contours for the to- 
tal dark matter (black) and the dark matter accreted from 
the four most massive satellites (black dotted), in MWl, 
II204, H258, and H258dark. All of the simulations that in- 
clude the baryons produce near-spherical, slightly oblate, 
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Figure 2. (a-c) The decomposed distribution of rotational velocities at the solar neighbourhood {7 < R < 8kpc; \z\ < 2.1kpc) for 
three simulated Milky Way mass galaxies MWl, H204 and H258. The lines show a double Gaussian fit to the dark matter (black; black 
dotted), the dark matter accreted from the four most massive disrupted satellites (green, blue, magenta, cyan), and the sum of all dark 
matter accreted from these satellites (black dotted histogram). The best fit double Gaussian parameters are marked in the top left, along 
with the redshift, z. (d) As (a-c), but for the galaxy H258 simulated with dark matter alone, (e-h) The decay in radius r as a function 
of redshift z of the four most massive disrupting satellites in MWl, H204, H258 and H258dark. Where less than four lines are shown, 
these satellites accreted at redshift z > 3. The dotted sections show the evolution of the most bound particle in the satellite after the 
satellite has disrupted, The decay in angle to the Milky Way stellar disc, 9 as a function of redshift z of the four most massive 

disrupting satellites in MWl, H204, H258 and H258dark. 



final dark matter density distributions that have symmetry 
axes parallel to the stellar disc. This agr ees well with pre- 
vious numerical results in th e literature l|Kazantzidis et ahl 
I2OO4I : iDebattista et all l2007l l. Indeed, observations of the 
Sagittarius stream of stars suggest that the Milky Way has a 
near-spherical halo with axis ratio a/c ^ 0.9 w ithin 50 kpc 
llJohnston et al. I I2OO5I : iFellhauer et al.l [2OO6I ). Of our three 
simulated galaxies, only the most flattened case, H204 with 
an axis ratio a/c ~ 0.77 within ~ 50 kpc, is inconsistent with 
the Milky Way. By contrast, without the baryons, the dark 
matter halo in H258dark is highly triaxial (Figure [3jh)). 
Indeed, the typical dark matter halo in ACDM - in simu- 
lations that model the dark matter alone - is triaxial, pro- 



late, and inconsiste nt with observations of the Milky Way 
l|Maccio et al.ll2007l ). 

As noted above, both H258 and H258dark have 
mergers at redshift 2 ~ 1. In H258 this merger produces 
a dark disc, but in H258dark it does not. Since mergers of 
this mass ratio define the post-merger plane of the disc, the 
difference cannot be due to disc plane dragging. Instead the 
difference is due to the halo shapes. In a static oblate po- 
tential like that in H258, particles conserve the z-component 
of their angular momentum vector. Any sense of rotation 
established in the merger will be preserved. By contrast, 
in a static tnaxial potential like that in H258dark, parti- 
cles moving on regular orbits do not explicitly conserve any 
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Figure 3. (a-d) The mean streaming velocities for the young stars (red) and gas (red dotted) in a slice \z\ < 1.1 kpc for MWl, H204, 
H258 and H258dark. Over-plotted are the circular velocity profiles (computed from the enclosed mass) for the dark matter (black), the 
sum of dark matter accreted from the four most massive satellites (black dotted), and the sum of dark matter accreted from the inner 
10% of the four most massive satellites (black dashed), (e-h) Projected density contours for the dark matter (black) and the sum of dark 
matter accreted from the four most massive satellites (black dotted), (i-1) Density of stars (red), stars+gas (red dotted), dark matter 
(black) and dark matter accreted from the four most massive satellites (black dotted and dashed as in (a-d)) as function of height z in 
a slice 7 < i? < Bkpc for MWl, H204, H258 and H258dark. For MWl, H204, H258, we overplot the o bserved star (solid blue c ircle) 
and star-l-gas (open blue circle) stellar density in the solar neighbourhood for the Milky Way (taken from lHolmberg fc Flvnrill200(]| ) . For 
H258dark, we show the profile for the alignment given in panel (h) (black), and one at 90 degrees to this (blue). 



component of their angular momentum vector. Their orbit 
planes precess and any sense of rotation established in the 
merger is rapidly lost. As a result, the final velocity distribu- 
tion at the solar neighbourhood in H258dark is a Gaussian. 
This demonstrates that in addition to disc plane dragging, 
the near-spherical halo that results once the baryons are in- 
cluded is a vital ingredient in the formation and survival of 
a dark matter disc. 



3.3 The mass and spatial distribution of the dark 
disc 

The precise properties of the dark disc depend on our choice 
of decomposition. Here, we define the dark disc as the mate- 
rial that originates from the innermost 10% of the four most 
massive merging satellites. This gives a reasonable match 
to the highly rotating Gaussian component in our double 
Gaussian fit to the local dark matter velocity distribution 
(see 

Figure [3j a-d) shows the mean streaming velocities {v^ 
averaged in bins AR — 1 kpc) for the stars (red) and gas 
(red dotted) in a slice \z\ < 1.1 kpc for MWl, H204, H258 
and H258dark. Over-plotted are the circular velocity pro- 
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Figure 4. (a-d) A comparison of the accreted star (red) and accreted dark matter (black) velocity distributions vji (left), (middle) 
and 1)2 (right) at the solar neighbourhood {7 < R < 8kpc; \z\ < 2.1kpc), for MWl, H204 and H258. We plot the sum of material accreted 
from the four most massive disrupted satellites. Also plotted is the dark matter coming from the inner 10% of these satellites (dotted 
lines). The dispersions are given in the top left corner of each plot in km/s. 



files (computed from the enclosed mass) for the dark matter 
(black), the sum of dark matter accreted from the four most 
massive satellites (black dotted), and the sum of dark mat- 
ter accreted from the inner 10% of the four most massive 
satellites (black dashed). Using our above decomposition, 
the dark disc contributes at most half of the rotation curve 
at the solar neighbourhood (in H258). Since the rotation 
curve goes as the square root of the mass, this is a quar- 
ter of the mass. The other extreme is MWl where the dark 
disc contributes a quarter of the rotation curve at the solar 
neighbourhood, which is 1/16 of the mass. 

Figure [3ji-l) shows the density of stars (red), dark mat- 
ter (black) , dark matter accreted from the four most massive 
satellites (black dotted) , and dark matter accreted from the 
inner 10% of the four most massive satellites (black dashed) 
as a function of height z in a slice 7 < R < 8kpc for MWl, 
H204, H258 and H258dark. In MWl, H204 and H258 the 
density profile is peaked, indicating a disc-like structure. In 



H258dark, the profile depends on our choice of alignment. In 
its eigenframe, p{z) is fiat, as expected for a more spheroidal 
dark matter distribution. In a frame rotated 90 degrees to 
this, however, the distribution is peaked similarly to H258 
(blue line; Figure |3][1)). However, there is no measurable ro- 
tation in the dark matter distribution for H258dark for 
any alignment. The peaky p{z) distribution in H258dark is 
the result of a fiattened halo, not a dark disc. For our above 
decomposition, the dark disc in MWl, H204 and H258 is a 
structure that rotates in the same sense as the stellar disc 
and is fiattened in the same plane as the stellar disc - hence 
the name 'dark disc'. 



3.4 Observing the dark disc 

Figure 3] shows a comparison of the star (red) and dark 
matter (black) velocity distributions at the solar neighbour- 
hood {7 < R < 8kpc; l^l < 2.1 kpc) of material accreted 
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from the four most massive disrupted satellites in MWl, 
H204 and H258. Also plotted is the dark matter coming 
from the inner 10% of these satellites (dotted line). The 
accreted stars comprise 5, 40 and 30% of the stars in the 
sola r neighbourhood fo r MWl, H204 and H258, respectively 
(see iBrook et al] l2009l for a study on the decomposition of 
these into bulge, disc and halo stars). The accreted stars 
are colder than the full accreted dark matter distributions, 
but closer to the inner 10% dark matter distribution. The 
best-fit parameters for the dark disc (see Figure [2ja-d)) lie 
somewhere in-between the stellar distributions and the in- 
ner 10% dark matter distributions. If future surveys of our 
Galaxy can disentangle accreted stars in the Milky Way 
thick disc from those that formed in-situ, then we will be 
able to infer, through numerical modelling, the velocity dis- 
tribution function of the dark disc from these stars. We will 
investigate this in more detail in future work. 



4 DISCUSSION: THE IMPORTANCE OF 
BARYONS 

The simulations MWl, H204, H258 and H258dark demon- 
strate the importance of including the baryons in simula- 
tions of the local dark matter distribution. Once the stars 
and gas are included, the dark matter halos become signif- 
icantly rounder, while the stellar and gas disc biases the 
accretion of massive satellites, dragging them towards the 
disc plane. Both effects aid the formation of a dark disc and 
lead to solar neighbourhood velocity distributions that are 
anisotropic, and better fit in by a double Gaussian than 
a single Gaussian. 

It is not clear which of MWl, H204 or H258 provides the 
best match to our Galaxy. The Milky Way has evidence for 
several mergers since z ~ 1. The Sagittarius dwarf recently 
fell i n on a polar orbit l|lbata et al.l 1 19941 : iMaiewski et al.l 
I2OO3I ) there is evidenc e for a merger in the plane of the 
disc ijConn et al.ll2007l ). and another ^ 5 Gyrs ago from 
the chemis try and kinematics of st ars in the solar neigh- 
bourhood (|Feltzing fc Bensbvl [ioosl ') . while the Milky Way 
thick disc itself is evidence for a ^1:10 merger at z ~ 1 
iRead et al.ll2008l : lKazantzidis et ai]|2007l '). MWl had no sig- 
nificant mergers after 2; ~ 2 and is likely over-quiescent as 
compared to our Galaxy. By contrast, both H204 and H258 
appear to have merger histories that are too extreme as com- 
pared with the Milky Way. This results in large fractions 
(30 — 40%) of accreted stars at the solar neighbourhood that 
are inconsistent with observations of the Milky Way. (As- 
suming all stars in the Milky Way thick disc a re accreted 
gives an upper bound of 12% for our Galaxy (jjuric et al.l 

mi)-) 

However, differences between the simulated galaxies and 
the Milky Way could also be a result of numerical limi- 
tations. All three state-of-the-art simulations have a force 
softening of ~ 300 pc, which is larger than th e measured 
scale height of the Milky Way thin stellar disc jjuric et al.l 
|2008| ): spurious angular momentum loss occurring in under- 
esolved progenitors may have still affected the final mass 
distribution even a t the high mass resolution employed here 
jMaver et al.ll200^ : the sub-grid physi cs scheme is a source 
of systematic error l|Maver et al.|[200^ : and the smoothed 
particle hydrodynamics technique employed does not cor- 



rectly resolve fluid instabilities and mixing l|Agertz et al.l 
l2007i i. 

Owing in part to these limitations in the numerics, and 
in part to the fact that the initial conditions were not chosen 
specifically to reproduce the structure and assembly history 
of the Milky Way suggested by observational constraints, the 
comparison with the Milky Way has to be done with caution. 
Indeed, the simulated galaxies more closely resemble L* Sa 
galaxies than a Milky Way- like Sb/Sc gal axy, as shown by 
their location in the TuUy Fisher relation l|Governato et all 
l2008f ). They have a large r bulge-to-disc ratio relative to the 
Milky Way (|Maver et al.l i,2008). ~ 0.5 -0.65 (measured 
both photometricaly and kin ematically, see lGovernato et al.l 
l2007l : ]Governato et al.ll2008l ') rather than the 0.2 - 0.3 mea- 
sured for the Milky Way, a more massive stellar halo and 
a lower gas fraction in the disc. As a result, while the total 
baryonic mass of these galaxies is comparable to that of the 
Milky Way, the stellar disc mass is a factor of 2-3 lower than 
that of the disc of the Milky Way, and the gas mass in the 
disc is lower by an even greater factor. The lower disc mass 
combined with the effect of softening, that tends to smear 
out physical density, yields a stellar volume density lower 
by a factor of ~ 3.5 compared to the stellar density of the 
Milly Way at the solar neighborhood and a total baryonic 
(star-|-gas) density in the disc plane lower b y a factor ~ 10 
at the same radius jHolmberg fc F Ivnn' 2000; and see Figure 
[3ji-l)). Due to the same structural differences, in particular 
to the larger bulge-to-disc ratio, the mean streaming velocity 
also falls with radius faster than in the Milky Way (Figure 
[Sja-d)). This will cause disc-plane dragging - and therefore 
the mass of the dark disc - to be underestimatec|f|. 

Despite these numerical and observational limitations, 
our qualitative results are robust. This is because dark disc 
formation is a gravitational process that requires just three 
key ingredients to be correctly modelled: (i) having a stel- 
lar/gas disc in place at high redshift as observed in real 
galaxies; (ii) having a dark matter halo that is oblate and 
aligned with the stellar disc; and (iii) having a merger his- 
tory given by our current ACDM cosmology. All three are 
satisfied by our simulations. 

With the above limitations in mind, we tentatively sug- 
gest that the Milky Way is intermediate between simulations 
MWl and II204. Future simulations that better resolve the 
Milky Way and its satellites and that are tuned to repro- 
duce the Milky Way's merger history will provide improved 
constraints. 



This is not likely to be the main source of error, however, since 
instabilities mostly affect the interstellar medium that is in any 
case poorly resolved. 

* There is a competing effect that has the opposite sign: the merg- 
ing sateUites are themselves over-concentrated due to the same 
numerical Umitations. This causes the dark matter, through adi- 
abatic contraction, to be over-concentrated also. This has two 
effects. Firstly the satellites are more resilient to tides and more 
able to get down into the disc. Secondly, there is more dark matter 
at small radii within the satellites. Both effects will act to overes- 
timate the dark disc mass. However, even without any baryons, 
H258dark has significant contributions from its accreted satellites 
at small radii (compare Figures [Sfc) and (d), dotted and dashed 
lines). This suggests that this effect is small. 
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5 CONCLUSIONS 

Predicting the local phase space density of dark matter is 
central to efltorts to directly detect dark matter, both to 
motivate detector design and to interpret any future signal. 
Previous efforts to estimate this have used simulations that 
model the dark matter alone. In this paper, we used three 
state-of-the-art ACDM cosmological hydrodynamic simula- 
tions of Milky Way mass galaxies to include the stars and 
gas self-consistently for the first time. Our main findings are 
as follows: 

(i) Once the stars and gas are included, the dark mat- 
ter halos become significantly rounder, while the stellar/gas 
disc biases the accretion of massive satellites, dragging them 
towards the disc plane. Both effects combined lead to the 
formation of a dark matter disc with 0.25 — 1.5 times 
the halo density at the solar neighbourhood, in excellent 
agreement wit h a previous estim ate obtained using different 
methodology l|Read et al.ll200§ '). The resulting solar neigh- 
bourhood dark matter velocity distributions are anisotropic, 
and better-fit in by a double Gaussian than a single Gaus- 
sian. The highly rotating second Gaussian component is the 
'dark disc'. We found a rotation lag compared to the Milky 
Way stellar disc in the range: wiag ~ — 150km/s. 

(ii) Our three Milky Way mass galaxies were chosen to 
span a range of interesting merger histories, not to be precise 
replicas of the Milky Way. We found that more late mergers 
gave rise to a more significant dark disc. The Milky Way 
likely has a dark disc intermediate between our most quies- 
cent galaxy, and our second most quiescent galaxy. This sug- 
gests a local dark disc density of pddisc/phalo ~ 0.23 — 1 
and rotation lag with respect to the Milky Way's stellar 
disc of wiag = — 150km/s. Increased resolution and cosmo- 
logical models that better capture the structural properties 
of the Milky Way will give improved constraints. For me- 
dian values of pddisc/phalo ~ 0.5 and uiag ~ 50km/s, 
the dark disc boost s WIMP capture in the Earth and Sun 
l|Bruch et al.1 l2008d : iBruch et aLlbOOSal ). enhances the an- 
nual modulation signal, and leads to distinct variations in 
the flux as a function of recoil energy that a llow the WIMP 
mass to be determined l|Bruch et al.ll2008bl ). 

(iii) An accreted disc of stars forms concurrent with the 
dark disc and shares similar kinematics. If future surveys 
of our Galaxy can disentangle accreted stars in the Milky 
Way thick disc from those that formed in-situ, then we will 
be able to infer, through numerical modelling, the velocity 
distribution function of the dark disc from these stars. We 
will investigate this further in future work. 

We conclude that it is vital to model the baryons when 
calculating the local phase space distribution of dark matter 
in the Milky Way. 
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